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Pressure effects on tfie cfiarge and spin dynamics in tfie biiayer manganite compounds 
La2-2xSri-|-2xMn207 are studied tfieoreticaify by taking into account tfie orbitaf degrees of free- 
dom. Tfie orbitaf degrees are active in tfie fayered crystaf structure, and appfied fiydrostatic 
pressure stabifizes tfie 3d^2_y2 orbitaf in comparison witfi 3do,22_^.2. Tfie cfiange of tfie orbitaf 
states weakens tfie interfayer cfiarge and spin coupfings, and suppresses tfie tfiree dimensionaf fer- 
romagnetic transition. Numericaf resufts, based on an effective ffamiftonian wfiicfi incfudes tfie 
energy fevef difference of tfie orbitafs, sfiow tfiat tfie appfied pressure controfs tfie dimensionafity 
of tfie spin and cfiarge dynamics tfirougfi cfianges of tfie orbitaf states. 
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Tlie colossal magnetoresistance (CMR) observed in 
several manganese oxides has recently attracted con- 
siderable attention from both academic and technolog- 
ical viewpoints. In order to improve the technical rele- 
vance, several material syntheses and processing meth- 
ods based on the manganites have been attempted. One 
of the convincing candidates for the large MR mate- 
rials is the manganese oxides with layered perovskite 
structure Lai_xSri_|_xMn04 and La2-2xSri_|_2xMn207. 
[1] The layered compounds consist of the single- or 
double-layered Mn02 planes and the intervening insu- 
lating planes stacked along the c-axis. Therefore, they 
are recognized as an intrinsic ferromagnetic metal (FM)- 
insulator (I)-FM multilayered system. 

The unique characteristic of the MR effect observed 
in Lai 4Sri 6Mn207 is interpreted as the tunneling mag- 
netoresistance. [2] The transport and magnetic proper- 
ties in this compound are characterized by two temper- 
atures, [2] T^'„^ and T^„^ (T^''^^ > T^^x), where the 
resistivities in the a&-plane (pah) and along the c-axis 
(pc), respectively, become maximum. In the region above 
"^maxj both pab and Pc increase with decreasing temper- 
ature. The remarkable feature appears in the interme- 
diate temperature range {T^^^ > T > T^^^), where 
Pah decreases with decreasing temperature, although pc 
still exhibits insulating behavior. A large MR value 
(pc(H = 0)//9c(i?=50kOe) ~ 10000%) is observed in this 
temperature range. It is interpreted that around T^^^ 
the ferromagnetic spin correlation in the a&-plane starts 
growing up. On the other hand, spins are almost ran- 
domly oriented along the c-axis . Below T^^^, the three- 
dimensional ferromagnetic correlation develops and the 
resistivities along both directions rapidly decrease. 

Recently, the pressure effects on the transport and 
magnetic properties in Lai 4Sri 6Mn207 are studied by 
Kimura et al. [3] It is shown that applied hydrostatic 
pressure up to llkbar suppresses the steep drops of p^ 
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and Pah around T^^ix and preserves the two-dimensional- 
like charge and spin dynamics down to low temperatures. 
In other words, the magnitudes of the interplane mag- 
netic and charge couplings are reduced by the applica- 
tion of pressure. As a result, the applied magnetic field 
promotes the transport along the c-axis and results in a 
larger MR effect than that obtained under the ambient 
pressure. This is highly in contrast to the manganites 
with three-dimensional perovskite structure, where the 
applied pressure increases the ferromagnetic transition 
temperature, and it is interpreted as the enhancement 
of the ferromagnetic double exchange interaction due to 
the pressure. [4] On the other hand, it is difficult to un- 
derstand the confinement-like behavior of the spin and 
charge dynamics observed in the layered compound due 
to the applied pressure from the conventional double ex- 
change scenario. 

In this letter, we study, theoretically, the pressure ef- 
fects of the transport and magnetic properties in layered 
perovskite manganites by taking into account the orbital 
degrees of freedom. It is well known that in manganites 
with three-dimensional perovskite structure the Cg elec- 
tron in the Mn^"'" ion has orbital degrees of freedom, in 
addition to charge and spin. The orbital plays several 
important roles in the spin and charge dynamics in man- 
ganites, as several authors have pointed out. [5-13] It is 
undoubted that the orbital degrees are also active in the 
layered manganites of the present interest. Actually, in 
recent X-ray diffraction experiments the orbital ordering 
has been directly observed in Lai sSro 5Mn04, [14] which 
supports the CE-type spin structure [15] known in the 
cubic compounds. [16,5] Because the orbital degrees of 
freedom strongly couples with the lattice, the structural 
change due to the applied pressure results in a modifi- 
cation of the orbital state. To confirm the change of the 
orbital state, we calculate the Madelung potential for the 
3dj,2_y2 and 3dg^2_^2 orbitals, and show that the 3dj,2_y2 
orbital is relatively stabilized with increasing pressure. 
By considering the energy-level difference between the 
two Cg orbitals, we derive the effective Hamiltonian de- 
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scribing the low-energy electronic states, and investigate 
the spin and orbital phase diagram. The numerical re- 
sults are satisfactory to explain the suppression of the 
three-dimensional ferromagnetic state experimentally re- 
ported. 

First, let us consider the orbital states under pressure. 
The crystal structure in Lai 2Sri 8Mn207 under the hy- 
drostatic pressure up to 6.3 kbar has been studied by 
the neutron diffraction experiment. [17] It is reported 
that the absolute value of the c-axis compressibility is 
almost three times larger than that in the a&-plane in 
the entire temperature range. In the MnOe octahedron, 
the shortening of the Mn-0 apical bond is remarkable 
(Ad/d ~ 1% for 6.3kbar) in comparison with that of 
the in-plane Mn-0 bond (~ 0.1%) in the ferromagnetic 
phase. The in-plane Mn-O-Mn bond angle remains un- 
changed (~ 180°), indicating no tilting of the octahe- 
dron. In order to estimate the orbital state under the 
above structural change, we calculate the Madelung po- 
tential which controls the energy levels of the Cg orbitals. 
The conventional Ewald method is adopted using the 
structural data at lOOK. The Madelung potential for the 
3dj;2_y2 and 3(i3^2_r2 orbitals is given by V(ro ± rdx) 
and (V(ro + rdz) + V(ro — rdz))/2, respectively. Here, 
T-Q is the position of the Mn ion, rd(= 0.42A) is the radius 
where the radial charge density of the Mn 3d orbital is 
maximum [18] and x and z are the unit vectors along the 
in-plane and inter-plane Mn-0 bonds, respectively. The 
potential difference between the 3dj;2_y2 and 3(i3^2_r2 or- 
bitals for the electron is obtained as (Vx2_y2 —Vsz2-r2) = 
0.0313eV (0.48kbar), 0.0046eV (3.1kbar) and -O.OOOTeV 
(6.3kbar). With increasing pressure, the potential for 
the — orbital is relatively decreased. Although the 
crystal structure was measured from 0.48kbar to 6.3kbar, 
we suppose that the relative energy stabilization for the 
x^ —y^ orbital amounts to an order of O.leV at ll.Okbar 
where the large pressure effect appears in the resistivity. 
The above changes of the potential are mainly attributed 
to the anisotropic change of the Mn-0 bond length in the 
octahedron. We confirm that the hydrostatic pressure 
acts like an uniaxial pressure applied along the c-axis in 
terms of the stabilization of the x^ — y^ orbital. 

Next, we study the spin and charge dynamics by tak- 
ing into account the change of the orbital state due to 
the applied pressure. We extend the model Hamiltonian, 
previously derived for the manganites with cubic struc- 
ture, [10] to the case of the layered compounds under 
pressure. For simplicity, we consider the uniaxial three- 
dimensional crystal structure consisting of Mn ions, in- 
stead of the actual layered structure. In this model, 
the bilayered Mn02 planes are represented by a single 
plane where the Mn ions form a square lattice. In each 
Mn site, we introduce two Cg orbitals and a localized 
spin (5'*=!') with S = 3/2 as the t2g spin. The differ- 
ence of the energy levels between a(= 3z^ — r^)- and 
h(= x^ — t/^)-orbitals is represented as A = — £j. The 
Hund coupling (Jh) between Cg and t2g spins and the 
electron-electron interactions (U,U',I) between Cg elec- 
trons are introduced. Furthermore, between the nearest 
neighboring i and j sites, the electron transfer (t]^ ) for 



the Cg electrons and the antiferromagnetic interaction for 
the t2g spins are considered. Among the above energy 
parameters, the electron-electron interactions U, U' give 
the largest energy scale. Therefore, we derive the effec- 
tive Hamiltonian to describe the low energy electronic 
state by the perturbational calculation with respect to 
(the transfer energyj^/(the electron-electron interaction) 
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The first term 



describes the Cg electron system and is given by 

i <z'j>(777' 

+ Hu + Hu'-i + Hu'+i , (1) 



where di^cr is the annihilation operator excluding the 
double occupancy, and Ti = \ Y^a-y-y' d\^a'^ii''^ii' a is the 
pseudo spin operator for the orbital. tJJ is represented 
a^^ttn;; where the numerical factor de- 
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pending on the orbital. [10] The energy-level difference 
between the two orbitals is represented by A Ti^ . The 
second term and the last three terms in eq.(l) are named 
as the 't-term' and the 'J-term', respectively, on the anal- 
ogy of the t-J model. Among the three 'J-terms', Hiji_j 
is the leading term and its explicit formula is given by 
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2 (t^'' +f,;;')(f.n,+n,+ + f+n,.n,.) 
+ 2 foit^f +t'if)(ni+nj. + ni.nj+) 

-t^"t^|4/o(T,+T,_+T,_T, + ) 

X ((/o + f-)nj+ - (fo + /+)ni-) 



X ((/o + f-)ni+ - (fo + /+)«»-) 



,(2) 



with Si = \ J2ai3-y A-ya^^apdi^fj and ni± 
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U'-I "^"^ j± — u'-i±A which are the inverse of 
the intermediate energy including the energy difference 
between a and h orbitals. The second and third terms in 



Hgf f are shown by 
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where J'^f is the antiferromagnetic superexchange inter- 
action between the nearest neighbor t2g spins. The uni- 
axial crystal structure is reflected in the anisotropy of the 

transfer intensity t]^ and the superexchange interaction 

t . . . 

JjJ", and their ratio is represented as i? = toi i^^ /ioi,i+z 

= ^{A%xlA%z)- For simplicity, t^i^i^x and J*^-!^^ are 
represented as to and J'^f^ respectively. The effects 
of the applied hydrostatic pressure is described by the 
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change of the energy level difference A. When we as- 
sume as A = and ij = 1, the Hamiltonian is reduced 
to the isotropic one. [10] In the following numerical calcu- 
lation, we use the model which is expanded with respect 
to A/[/ up to the first order. Because the A depen- 
dent part of the 'J-term' shows 0(At^/U^), the main 
contributions of the level difference comes from the first 
term in eq. (1). When the transfer matrices are as- 
sumed to be tfj = Sabtoij in eq. (1) at a; = 0, we ob- 
tain the effective Hamiltonian, where the orbital space 
is isotropic, [19,20] plus A^^. Tj,j. In this model, the 
ground state changes from the well-known ferromagnetic 
spin ordering (spin:F) accompanied with the antiferro- 
magnetic orbital ordering (orbitahAF) to the spin:AF 
with the orbitahF at A ~ 3I(t^/U'y. In other words, 
the spin ferromagnetic state becomes unstable by lifting 
the orbital degeneracy. 

Based on the above Hamiltonian, we calculate the 
spin and orbital phase diagram in the mean field ap- 
proximation. We consider four types of the spin and 
orbital orderings in the cubic eight Mn unit; ferro- 
magnetic (F) ordering and layer- type (A), rod-type (C) 
and NaCl-type (G) antifferomagnetic(AF) orderings. In 
the 'J-term', (Siz) and (rij) are adopted as the or- 
der parameters. Because of the anisotropy in the or- 
bital space, we introduce the rotating frame in the 'J- 
term' where the orbital state is described by the an- 
gles (6''-^\6i^^) for each orbital sublattice as \0^*^) = 
cos(6'fV2)|3z2 - r2) - sm(0^p /2)\x^ - y^). Then, we 
adopt (Tj'^)(= cos 6'j-*-'(Tj'^) + sin O'f ^Tij;)) as the order 
parameter for the orbital. By the same token, the rotat- 
ing frame is adopted in the 't-term' by introducing the 
unitary matrices in the spin and orbital spaces, dia-y is 
transformed as dicr-y = h\z\^Jzf^ where hi is a fermion 

operator describing the hole motion, and zl^J and zf^ 
are the elements of the unitary matrix in the spin and 
orbital frame, respectively. Then, 't-term' is rewritten as 

-E<,,>/*I/*.(E.4:^*-?;)(E,y4?:^?,^'4^0, [11]- In 
order to qualitatively discuss the doping dependence of 

the phase diagram in the low doping region, we adopt 
the mean field approximation and replace {h\hj) by the 
carrier number x. Iz"/}* z^f}) and (z.-*-'*t77 z"i^}i) are rep- 
resented by the angles in the spin and orbital spaces, 
respectively. When we neglect the azimuthal angle in 
the spin space, '^^ri^ia^* ^jJ) reproduces the form as 
cos((6'j-^-' — 6'j^-')/2) derived by Anderson and Hasegawa. 
[21] The spin and orbital phase diagram at T = as 
a function of x and .P^^ is obtained by minimizing the 
energy. 

In Fig. 1, we present the numerical results of the spin 
and orbital phase diagram. In the numerical calcula- 
tion, we only consider the leading term Hjj'-i among 
the three 'J-terms'. The parameter values are chosen 
as U' = 3.0^0, I = l.Oto and R = 2.0. First, we focus 
on the A = 0.0 case (Fig. 1(a)). It is noted that in 
the case of i? = 1, the present calculation well repro- 
duces the results obtained by the exact diagonalization 
method based on the effective Hamiltonian at a; = [12] 
and the Hartree-Fock calculation in the region of finite 
X. [13] Therefore, the lack of the spin:C phase x = 0.0 



between the spin:A and spin:G phases is attributed to 
the anisotropy in tJJ and J*J' . With increasing the 
carrier concentration from x = 0.0, the region of the 
ferromagnetic phase gradually grows up, and becomes 
maximum around x = 0.1. With further increasing the 
concentration, this phase disappears aX x = 0.2 and the 
A-type AF phase becomes stable in the large parameter 
region. The stabilization of the A-type AF in the region 
of higher x originates from the x'^ — orbital ordering 
which is most favorable to the gain of the kinetic energy 
of the doped carriers, as indicated by the Hartree-Fock 
calculation. [13] Although the superexchange and dou- 
ble exchange interactions contribute to the ferromagnetic 
phase which appears around 0.0 < x < 0.2, the former is 
mainly dominant. The orbital state in the ferromagnetic 
phase is gradually changed from [o'f ^ = 90° , ^ = 270°) 
at around x = 0.0 to (6*^*^ = 120°, 6*^*^ = 240°) around 
X = 0.2. When A is introduced (Fig. 1(b)), re- 
markable changes appear around the boundary between 
the spin:F and spin:A phases, that is, the large region 
of the ferromagnetic phase is replaced by the spin: A 
phase. Around the boundary, spin:F accompanied with 
[e^P = 120°, 6*^*^ = 240°) energetically competes with 
spin:A with the x^ — j/^-type orbital ordering. There- 
fore, the small value of A relatively stabilizes the spin:A 
phase. The suppression of the three-dimensional ferro- 
magnetic phase originates primarily due to the following 
two reasons: i) the magnitude of the ferromagnetic su- 
perexchange interaction is reduced by lifting the orbital 
degeneracy, and ii) both super and double exchange in- 
teractions in the c-direction are prevented in the x^ — 
orbital ordering. On the other hand, the spin:F phase 
around a; = 0.0 is still robust in the case of A = O.lto, 
because this phase does not adjoin the other spin phase 
accompanied with the x^ — j/^-type ordering. In Fig. 2, 
we present the A dependence of the angles (6'^*'', 6*2*'') 
in the orbital space at X = 0.16 and = O.OlT^o- 

As the magnitude of A is increased, the pseudo-spins 
gradually cant from (6*^*^ = 120°, 6*^*^ = 240°) toward 
(6*^*) = = 180°) (a;2 - j/2 orbital) within the spin:F 
phase, and finally the spin:A appears at A = 0.075to- It 
is clearly shown that the applied pressure tends to align 
the orbital pseudo-spin ferromagnetically, like an applied 
magnetic field on the spin space. 

It is shown by the numerical calculations that the 
spin:F phase near the boundary is replaced by the spin:A 
phase by introducing A. In this phase, the electron 
transfer is prohibited in the c-direction, and the inter- 
layer spin coupling is dominated only by the weak anti- 
ferromagnetic interaction jj^+z derived by the electron 
transfer in the Mn 3d-0 2p tt bond. As a result, both spin 
and charge couplings are considerably weak along the c- 
axis in comparison with those in the a&-plane. Also in 
the state where the ferromagnetic phase remains under 
the applied pressure, the weight of the x^ — y^ orbital 
becomes dominant as shown in Fig. 2, and the inter- 
layer coupling is weakened. Therefore, the above states 
are regarded as the two-dimensional ferromagnetic state 
with the weak interplane spin coupling, corresponding 
to the state observed in Lai 4Sri 5Mn207 under pres- 
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sure. Although our calculation is limited in the mean 
field phase diagram, where the qualitative discussion is 
limited in the low concentration region and the spin and 
orbital fiuctuations are neglected, the obtained result 
is sufficient to explain the confinement of the spin and 
charge dynamics within the a&-plane by the application 
of pressure. In our calculation, the actual layered struc- 
ture is replaced by the anisotropic three-dimensional one, 
as explained above. However, the ferromagnetic cou- 
pling within the bilayers is expected to remain under 
the applied pressure as the experimental results show; 
the length of the Mn-O-Mn bond linkage between two 
Mn02 planes shows weaker changes in comparison with 
that of the Mn-0 apical bond. Also, the present results 
are consistent with the temperature dependence of crys- 
tal structure in Lai 2Sri 8Mn207, where the Mn-0 apical 
bond length is sharply increased below the ferromagnetic 
transition temperature. [22] The applied pressure sup- 
presses an elongation of the apical bond and stabilizes 
the x'^ — orbital. As a result, the ferromagnetic tran- 
sition is restrained. It is noted that by application of 
uniaxial pressure along the c-axis, the two-dimensional 
properties are expected to be more significant in com- 
parison with those in the case of hydrostatic pressure. 

In summary, we investigate the pressure effects on the 
spin and orbital states in the double layered perovskite 
manganites. Because of the anisotropic change of the 
Mn-0 bond length, the change of the Madelung poten- 
tial by the application of pressure stabilizes the 3dj;2_y2 
orbital in comparison with the 3dg^2_r2 orbital. The 
spin and charge dynamics are restricted within the aft- 
plane under the 3dj;2_y2 orbital state. As a result, the 
three-dimensional ferromagnetic state is replaced by the 
two-dimensional one accompanied with the weak inter- 
plane spin coupling. The applied pressure in the layered 
compound is recognized as a tool to control the dimen- 
sionality of the spin and charge dynamics through the 
changes of the orbital. 
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Figure captions 

Fig. 1: The spin phase diagram with and without the 
energy level difference A between two Cg orbitals. The 
parameter values are chosen as U' = 3.0to, I = 1-Oto and 
R = 2.0. The dotted line in Fig. 1(b) shows the phase 
boundary between the spin:F phase and the other spin 
phases for A = 0.0 case. 

Fig. 2: The angles (6'^*-', 6'2*'')in the orbital space as 
a function of A. The parameter values are chosen as 
U' = 3.0to,I = l.Oto, = 0.017^0, X = 0.16 and 

R = 2.0. The angles 120°, 180° and 240° correspond to 
the dgj,2_^2, dj.2_y2 and dgy2_^2 orbitals, respectively. 



